Controlled engineering of extended states in disordered systems 
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We describe how to engineer wavefunction derealization in disordered systems modelled by tight- 
binding Hamiltonians in d > 1 dimensions. We show analytically that a simple product structure 
for the random onsite potential energies, together with suitably chosen hopping strengths, allows a 
resonant scattering process leading to ballistic transport along one direction, and a controlled coexis- 
tence of extended Bloch states and anisotropically localized states in the spectrum. We demonstrate 
that these features persist in the thermodynamic limit for a continuous range of the system param- 
eters. Numerical results support these findings and highlight the robustness of the extended regime 
with respect to deviations from the exact resonance condition for finite systems. The localization 
and transport properties of the system can be engineered almost at will and independently in each 
direction. This study gives rise to the possibility of designing disordered potentials that work as 
switching devices and band-pass filters for quantum waves, such as matter waves in optical lattices. 

PACS numbers: 71.30.+h, 72.15.Rn, 03.75.-b 



I. INTRODUCTION 

The electronic properties of most materials are deter- 
mined by their crystal structure or lack thereof. For 
purely crystalline materials, well established methods of 
condensed matter physics 1 - allow the nearly complete ex- 
perimental characterization, as well as theoretical de- 
scription, of the resulting electronic bands and associ- 
ated density of states (DOS), all the way to transport 
and thermal properties. In fact, our present under- 
standing of the underlying mechanisms allows the man- 
ufacture of tailored artificial crystal structures such as 
photonic £^ phononic^^ polaritonio&i or plasmonio^ 
lattices. While in the former two systems classical coun- 
terparts of electronic bands and transport properties are 
manifestly observed, in the latter two, electronic quasi- 
particle scattering is shown to lead to the formation of 
controllably engineered excitation bands and, in particu- 
lar, the gaps between these as is required for a multitude 
of applications ji ^— 

In a strongly disordered system, Anderson 
localization 1 ^ suppresses transport even in regions 
with a finite DOSi 14 ' 15 Particularly in low-dimensional 
systems, the so-called scaling hypothesis^ establishes 
the expectation that all states remain localized for non- 
interacting quasi-particles, and hence there seems little 
room for a similarly controlled "engineering" of bands 
of extended states. Nevertheless, such a complementary 
approach has already enjoyed some successes. Local 
positional correlation in a disordered material has been 
shown to lead to resonant scattering events generating 
extended states at isolated energies in the spectrum^ - — 
With much longer-ranged correlated disorder, even in 
one-dimensional (ID) systems, effective metal-insulator 
transitions can be inducedi^Ir— Similarly, certain dis- 



ordered configurations can lead to an optimization of 
the quantum interference process mediating excitonic 
transport in molecular networks ! 27 ' 28 In fact, a con- 
trolled disorder can induce highly selective transport 



properties 



29-32 



giving rise to materials with interesting 



new functionalities, as has been recently explored with 
photonic crystals M^— 

In this paper we show how to open a channel of per- 
fect transmission in an otherwise disordered system. Our 
approach uses a simple product structure for the random 
potential which, except for certain resonance conditions, 
leads to anisotropically localized states. We believe our 
results to be of particular relevance for ultracold atoms or 
Bose-Einstein condensates in optical lattices as these are 
ideal for studying disorder effectSi^Zr— In fact, very re- 
cently the first direct observations of localization of mat- 
ter waves in 3D disordered optical potentials have been 
reported ! 40 ' 41 Here, by engineering the underlying disor- 
der, our study shows how extended matter waves emerge 
from a background of localized states (cp. Fig. [T]). Fur- 
thermore, our disorder structure allows for independent 
tailoring of the transport properties of the system in each 
direction, and gives rise to energy coexistence of extended 
and localized states, which can be manipulated with a 
high degree of control. We prove our findings analyt- 
ically and we corroborate them by extensive numerical 
studies. 

In Section [TTl we describe the chosen xy-disordcr, re- 
stricting ourselves to 2D for presentational simplicity. We 
then deduce the form of the resulting DOS and show 
that the eigenstates exhibit localization. In Section ITLTl 
we discuss the conditions for the existence of extended 
states and show the appearance of a perfectly transmit- 
ting channel. Section IIVI presents a discussion on the 
engineering of the transport properties and spectral re- 
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FIG. 1: (Color online) Probability density IV'x.j/l °f eigenstates of an xy-disordered system with (a) extended in x and (b) 
localized behaviour. The parameters characterizing the potential in case (a) are m a = 0, W a = t, mp = 1, W 7 ^ = 1, 
i v = £ = 0.55488 and eigenenergy S/t = 1.7862. For (b) we have m a = 0, W a = t, m = 4, Wp = 1, m £ = 1, W 5 = 0.5 and 
eigenenergy E/t = —0.10995. In both cases the system length and width are L x = 100 and L y = 50 respectively. The left and 
bottom panels show respectively the distributions \xy\ 2 and |<?M 2 that construct the eigenstate according to (|13|) . The color is 
determined by the value of — log^ IV^.i/l 2 where (L) = ^/L X L V (in this scale, 2 is the average value of an extended state). A 
second order interpolation of the distributions is used to smooth the visualization. 



gions with mixed localized and extended states. In Sec- 
tion H we draw our conclusions. Details on numerical 
techniques and some lengthy derivations are provided in 
several Appendices. 

II. LOCALIZATION AND SPECTRAL 
PROPERTIES IN xy-DISORDER 

A. The model 

We work with a 2D tight-binding model described by 
the Hamiltonian 

L x L x 1 

U = Y1 c Uxc x + ( c i tc *+i + 4+i tc x), (!) 

x—1 x — 1 

on a lattice with L x x L y sites. Here, e x denotes the 
L y x L y Hamiltonian matrix acting in the y direction for 
each vertical arm at position x such that 

Ux,\ lx,i ••• ••• \ 

7z,l £x,2 lx,2 
e = 7^,2 £x,3 7*, 3 ^ 

: '■• '■• '•• 

: 7x,L„-2 Cx,L v -l Jx,L v -l 
V Jx,L v -l Cx,L y J 

and 4 = (4,i. 4,2. ■ • ■ . 4,l J . with c *,y (4, y ) the usual 
annihilation (creation) operators at the site with coordi- 
nates {x,y}. Also, t = 41 is the hopping along the x 



direction. The set {e x ,y} gives the on-site energies and 
lx,y is the hopping term in the y-direction between sites 
(x, y) and (x, y + 1). A pictorial representation of the 
lattice is shown in Fig.[2][a). 

The energies and states of T-L are the solutions of the 
Schrodingcr equation 

(£1 =*(¥*+! + **_!), (3) 

where ^ x = {ip x ,i,ip x ,2, ■ ■ ■ , ipx,L y ) contains the wave- 
function amplitudes in the x-th vertical arm of the system 
and E is the energy. Here and in the following, we shall 
assume hard wall (fixed) boundary conditions in x and 
y direction, although the formalism does of course work 
with periodic boundaries as well. 

B. The disorder 

We now introduce disorder into the system in the fol- 
lowing way. Let {<>,}, i / and {P y } y =i,...,L y be ran- 
dom, uncorrelated sequences with corresponding proba- 
bility distributions V(a) and V(/3) which for simplicity 
will be taken as box-distributions of widths W a , Wp , and 
mean values (a) = m a , (/3) = nip. The on-site random 
energies will be constructed from the product of these two 
sequences, i.e. 

^x,y — Q-xfiy- (4) 

Additionally, we also choose the hopping strengths in y 
direction to be random and determined by 

lx,y=U x £, y , (5) 
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FIG. 2: (Color online) 2D lattice with sy-disorder: (a) microscopic link-node representation. Different colors highlight different 
on-site energies at every site. The width of the links is proportional to the magnitude of the hopping terms. A twist in the 
link indicates a negative hopping term. The horizontal (a;) and vertical (y) diagrams show the sequences {a x } and {/3 y }, {£ y } 
respectively. Density plots for the on-site energies e x , y — a x /3 y and the vertical hopping terms "j x , y — for a 50 x 50 lattice 
with m a /t = ma = = and W a /t = Wp = W$ = 4 are shown in (b) and (c) respectively. A second-order interpolation of 
the random sequences is used in (b) and (c) to smooth the visualization. 



where {^ y } y= i l„-i is another independent random se- 
quence. For simplicity, we choose the elements {a x } to 
have dimensions of energy — measured in units of the 
hopping t — , thus {P y } and will be dimension- 

less. This choice of e XjV and j Xty above gives rise to 
patterns for the disordered energy landscape and the ran- 
dom (vertical) y couplings with a characteristic correla- 
tion in the x and y directions, as can be seen in Figs.[2jb) 
and[2Ic), hence the name 'ccy-disorder'. Let us remark 
that these choices, particularly for L y <C L x , resemble 
and generalise quasi-lD "ladder" models currently dis- 
cussed in the literature to describe electronic transport 
in DN A 42 i 43 and mcsoscopic devices 4^— However we 
believe that a faithful realization of our model can be 
implemented using optical potentials and matter waves, 
where single-site resolution and control has already been 
demonstrated 



C. Reduction to decoupled channels 



where P does not depend on the x coordinate. The ma- 
trix P can be diagonaliscd via p = U _1 PU, where U is 
the matrix whose columns are the orthonormal eigenvec- 
tors of P and p contains the eigenvalues pi, . . . ,Pl v in 
its diagonal. Performing the change to a new basis 



(7) 



where <& x = (j) x l \ <f> x 2 \ ■ . ■ , 4> x Ly ^ , we can reduce 
Eq. © to 



(El - a x p)$ x = t(* x+1 + * x _i) 



(8) 



This is simply a set of L y decoupled Schrodinger equa- 
tions, 



Equations Q and allow us to factorize the e x ma- 
trix as 



(E - a x p c 



(9) 
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(E - a xPLy )^ =t($& + <P ( tl 



each of which corresponds to a ID disordered channel 
with random on-site energies e x = a x p c for the c- 
th channel. The disorder is uncorrelated within each 
channel and the distribution is V{e) = V(a)/p c , which 
is again a box-distribution with mean p c m a and width 
IPcIW-q, for the c-th case. 



4 



D. Density of states 



From Eqs. ([9]), it follows that the energy spectrum of 
the 2D L x x Ly system can be obtained from the union 
of the different spectra for the corresponding decoupled 
ID channels. The spectrum will be determined by the 
properties of the random distributions {a^}, {P y }, and 
{£ y }- For finite L y the DOS per site of the 2D system 
can be written as 



g iD (E;m a ,W a ,mp,Wf,,{£y}) 



^Y,9( E ^m a APc\W a ), (10) 



where g(E;p c m a , \p c \W a ) is the DOS per site of the c- 
th decoupled channel, characterized by a disorder distri- 
bution of mean p c m a and width |p c |W a . In the limit 
L y — > oo, Eq. (fTU)) becomes 

g w {E-m a ,W a ,mp,Wp,{^ y }) = 



g(p;mp,Wp,{£ y }) g(E;pm a , \p\W a )dp, (11) 



where Q(p;mp,Wp,{£ y }) corresponds to the DOS per 
site for the eigenvalues of P. Due to the nature of xy- 
disordcr, the DOS of the 2D system follows from a con- 
volution of ID distributions. 

Eqs. ([TU|) and (TTT]) can be used to obtain the DOS 
numerically. Furthermore, since they only require cal- 
culations of 1-D distributions, they can be very effi- 
ciently combined with the functional equation formalism 
(FEF) to obtain the DOS of the system for L x = oo 
and either finite L y or L y = oo. A brief overview of the 
FEF is given in Appendix IA"1 Some examples of the DOS 
for cry-disorder are shown in Fig. [31 where we see that 
the global distribution of states is a superposition of the 
DOS from the individual channels. 




FIG. 3: (Color online) DOS g 2I) (E/t) and maximum localiza- 
tion lengths X(E) in the x and y directions for a;i/-systems as 
computed from the FEF. For panel (a) L x — oo and L y = 4 
for m a = t, W a = 2t. The dashed lines show g(E/t)/L y for 
the L y decoupled infinite ID channels that get added accord- 
ing to Eq. pop to produce g 2D (E/t). For comparison, the 
shaded histogram has been obtained from exact diagonaliza- 
tion of a (4 x 1000)-site system, averaging over 25 realizations 
of the a x sequence. The inset shows the fixed finite sequences 
{/3 y } and {£, y }. (b) L x = L y = oo with parameters mp = 4, 
Wp = 2, W a = 3t and m a = (solid), m a — At (dashed). In 
this case f s = £ = 1. The DOS was calculated using Eq. 
where the integral over p £ [1,7] was discretized using 1000 
points. Notice the different scale (right) for the localization 
lengths, obtained according to Eqs. (|14[) and (|15[) , 



E. Localization of eigenstates 

The states of the 2D system can be obtained from the 
eigenstates of the ID decoupled channels by undoing the 
change of basis Q, ty x = XJ& X . If we consider the n- 
th eigenenergy of the c-th decoupled channel with corre- 
sponding eigenstate {(f>lc }x=i,... ,L X > the vector & x will 
only have one non-zero component at the c-th position. 
It then follows that 

= u- (o,...,o,4 c,n) ,o,...,o) T 

= U (c U { x ' n \ (12) 

where U^ c ^ means the c-th column of U, which corre- 
sponds to the c-th eigenvector of the matrix P which we 
denote by {Xy C ^}y=i L y - The 2D eigenstates of ([TJ arc 



thus obtained as 

^y = Xy C) ^ n \ X Z\-- L * (13) 

y — 1, . . . , ljy 

where we must consider all ID eigenstates n = 1, . . . , L x 
for all the decoupled channels c = 1, . . . ,L y , giving the 
complete basis of L x x L y states in the original coordi- 
nates. Notice that all eigenstates of a certain decoupled 
channel c get multiplied by the same eigenvector of P. 

From Eqs. ([6]) and ([8]), and making use of known re- 
sults on Anderson localization in 1D^ we conclude that, 
for generic cases, the eigenstates (|13[) will be anisotrop- 
ically localized for all energies, with different localiza- 
tion lengths in the x and y directions. Localization per- 
sists although our model (|T|) only contains L x + L y — 1 
independent random on-site energies {e Xl i, ei. y , for x = 
1, . . . , L x , y = 2, . . . , L y }, as compared to the standard 
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Anderson model with L x x L y uncorrected values for 
e xy . An example of a localized eigenstate is shown in 

Fig.mb). 

For a given energy E, the transport properties of the 
system in the x direction, as L x — > oo, will be determined 
by the largest possible value of the localization length in 
this direction, A^ (E) . It then follows that 

X x {E)=max{\^(E)) , (14) 

^ J c— l,...,L y 

where the localization length for each channel c at en- 
ergy E is defined as the inverse of the corresponding Lya- 
punov exponent, X x c \e) = 1/t} x c \e)^ Additionally, lo- 
calization in the y direction will be determined by the 
Lyapunov exponent of the ID system defined by the P 
matrix only, \ v (p) = l/i} y (p) 1 for large enough L y . The 
localization length in y will change with the eigenvalue p, 
and thus all wavefunctions (|13|) obtained from the same 
decoupled channel [same p c in Eqs. ([5])] will exhibit the 
same y-sprcading. The relevant localization length in y 
as a function of the energy can be defined as 

\ y (E) = max{X y (p c )} pc/EeSc , (15) 




10 20 
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FIG. 4: (Color online) Full spectrum of Lyapunov exponents 
r]x(E) from TMM (dashed lines) and FEF (solid lines) for a 
system with a;?/-disorder (L x = oo, L y — 10) and parameters 
m a — At, W a = 2t, and £ w = £ = 1. The inset shows the 
10-element sequence considered for j3 y . The thick black line 
highlights the trajectory of the minimum Lyapunov exponent 
for each energy, i.e. the inverse of the maximum localization 
length, as computed by the FEF. The blue circles show every 
50th data point of the corresponding TMM result (errors are 
within symbol size). 



where 



III. DELOCALIZATION IN zy-DISORDER 



S r .= 



—— - 2t,p c m a + a—- 



21 

(16) 

denotes the boundaries of the energy spectrum of the 
c-th channel as L x — > oo, and a = sign(p c ). That is, 
at energy E, the relevant localization length in y will be 
the maximum of X y (p c ) over all channels whose spectrum 
includes E. 

The Lyapunov exponents for the decoupled ID chan- 
nels, in the thermodynamic limit, can be obtained numer- 
ically using the FEF (see Appendix IA"]l . which in combi- 
nation with (fT4"|) and (|15[) permits the calculation of the 
relevant localization lengths for a 2D system with generic 
icy-disorder, like those shown in Fig. [3] We emphasize 
that, irrespective of the detailed energy dependence of 
these localization lengths, all A values are finite and hence 
correspond to localized states. 

In order to confirm the validity of the calculation of 
the localization lengths, we carried out transport simula- 
tions in the x direction using the transfer-matrix method 
(TMM)±£ (see Appendix [Bj on Eq. © without perform- 
ing any decoupling. This technique is capable of giving 
the whole spectrum of characteristic Lyapunov exponents 
and thus automatically provides the largest decay (local- 
ization) length for an initial excitation travelling through 
the system. In Fig. 2]we show the whole Lyapunov spec- 
trum of an xy-disordered system (L x — > oo and finite 
Ly) obtained from TMM, and compare it with the results 
from the FEF on the decoupled channels. The excellent 
agreement observed confirms the correctness of the de- 
coupling transformation and of Eq. (|14[) . 



Despite the localized nature of eigenstates in generic 
xy-disorder, a genuine delocalization in the x-direction 
can be achieved as well. This is demonstrated in Fig. [5J 
where we show results from TMM calculations of the lo- 
calization length X X (E) for a system with L y = 10 and 
different values of £ y = £ (here we assume a constant £ for 
presentational simplicity only) . We see that the enhance- 
ment of the localization length is always found inside the 
region \E\ ^ 2. In fact, as shown in Fig. [51 for certain 
£ the localization length becomes comparable to L x in- 
dependently of the longitudinal disorder W a , signalling 
the appearance of extended states. This is in contrast to 
the W~ 2 dependence expected for quasi- ID systcmS ) 56 i 57 
which is found at most other £ values. 



A. Zero eigenvalue condition 

The reason for the existence of this band of extended 
states can be understood by returning to Eq. We 
see that if at least one of the eigenvalues {p c } of P is 
zero, then the corresponding ID equation in ([9]) no longer 
describes localization but rather an extended channel. As 
cachp c value depends on the distributions {(3 y } and {£j,}, 
the condition p c = can be controlled by fine-tuning 
these parameters in our 2D disordered system, therefore 
inducing a perfectly transmitting channel. 

Indeed, the relation between the {£ y } and {(3 y } val- 
ues can be derived straightforwardly. The characteristic 
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FIG. 5: (Color online) Energy dependence of the normalized 
localization length in the x direction obtained from TMM for 
a system of L y = 10, m a = and W a = 2t for different values 
of £y = £. For clarity only every 50th data point has been 
indicated by a symbol. Errors are within symbol size. For the 
chosen realisation of j3 y values (shown in inset), £ = 0.26042 
corresponds to p c = for one of the 10 channels. Note that 
the TMM calculations were terminated at a fixed value of 
L x = 10 6 and hence the X x values for £ = 0.26042 in the 
energy range \E\ ^ 2 do not diverge. The right frame shows 
the dependence of \ X /L X versus W a at energy E = 0. The 
green lines represent W^ 2 power laws. 



polynomial of P is 



trP 



• + f(P v , e,)p+(-l) iH detP = 0, 

(17) 

and if det P(/3 y , £ y ) = then this will correspond to at 
least one eigenvalue of P being zero. Assuming that the 
(5 V disorders have been fixed, the condition can be sat- 
isfied by choosing the remaining £ y values appropriately. 
For example, for L y = 2, 3, 4, the condition p c = reads 



Pi +fo 



PiMl 



/Wi 2 



(18) 
(19) 

/3i/3 2 /3 3 /3 4 , (20) 



respectively. In the most general case one could select 
arbitrarily L y — 2 values for £ y and then choose the re- 
maining one such that det P — 0. From here on, we shall 
restrict ourselves to the case where £ y = £ constant for all 
y. This simplification makes the hoppings in the y direc- 
tion constant within each vertical arm, i.e. 7. 



tC it 

should be clear, however, that this condition is in general 
not necessary for the existence of extended states. 

Because of the tridiagonal nature of P, its determinant 
involves only even powers of £, and it is a polynomial of 



order L y [{L y — 1)] for even [odd] 



L y . We can then con- 



sider only £ > without loss of generality. Therefore, the 
spectrum of £ satisfying detP(£) = contains at most 
Ly/2 or (L y — l)/2 independent positive real values^ For 
any of these resonant £ values^ (p c = 0), the spectrum of 
the otherwise disordered system has at least one perfectly 



conducting channel populated by L x extended states in 
the x direction, in the energy interval E € [— 2t, 2i\. 



B. Spectrum and stability of resonant channel 

The eigenstates of the resonant channel [with p c — 
1 Eq. flSJl] are Bloch states^ 



i)|2 _ 



L, 



1 



■ sm 



1 



(21) 



with energy E^ = 2t cos [nn/(L x + 1)], n = 1, . . ., L x , 

for hard-wall boundary conditions 4>^' n ^ = ^fc™! = 0. 
However, the 2D eigenstates given by Eq. (fl"5|) . are in 
general still localized in the y direction. An example 
of such states is shown in Fig. HJa). The presence of 
these extended states changes also the properties of the 
DOS of the system. According to Eq. ([TTJ]). the DOS now 
contains the contribution from the spectrum of a periodic 
chain, g{E) = Tr~ 1 {4t 2 ~E 2 )- 1 / 2 for E e [—2*, 2t], and its 
characteristic Van-Hove singularities at the band edges, 
as shown in Fig. EJa). 

In general, there is a clear correlation between the de- 
pendences of detP and \ x on £, as shown in Fig. [7] 
Therefore very large localization lengths occur in the 
vicinity of the resonant £ values, which will lead to ef- 
fectively extended channels in finite systems, even if the 
resonant condition is not precisely met. For the case 
considered in Fig. [71 corresponding to a quasi-lD system 
with L y = 10, we see that deviations from the resonant £ 
value up to ±10% still ensure localization lengths in the 
x direction of at least several thousands of lattice sites. 
More interestingly, we have also checked that 5% of ran- 
dom spatial deviations in each j x , y value — i.e. when 
the factorization §5§ for the vertical hoppings is weakly 
broken — , do have a similar effect only. This tolerance of 
the derealization channel is very significant, and makes 
our results relevant for potential experimental realiza- 
tions, where deviations from the theoretically obtained 
parameters are to be expected. Furthermore, the situ- 
ation becomes even more robust when the width of the 
system is increased, as we discuss below. 



C. The resonant condition for large L y 

The resonant condition requires p = to be an eigen- 
value of the matrix P, defined in Eq. For short L y 
this implies tuning the vertical hopping elements via £, 
for which several isolated values lead to the appearance 
of a transmitting channel. For large L y , however, the res- 
onant condition translates into whether the DOS for the 
system described by P satisfies g(p = 0) 7^ 0. If the latter 
is true, then p = is an eigenvalue as L y — > 00, and at 
least one of the decoupled channels of the system will not 
correspond to a disordered chain. Therefore the infinite 
system (L v — 00, L x = 00) will always exhibit a perfect 
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FIG. 6: (Color online) DOS and inverse localization length 
in the x direction for a;j/-disorder with resonant channel, ob- 
tained from FEF. (a) L x = oo and L y = 8 for m a = 2t, 
W a = 3t and £ = 0.44380. The dashed lines show g(E/t)/L y 
for the decoupled infinite ID channels, that get added accord- 
ing to Eq. I|10p to produce g 2D (E/t). The resonant channel 
contribution is highlighted in gray. The inset shows the fixed 
finite sequence {fi y }. (b) L y — L x = oo with parameters 
mp = 2, Wg — 2, W a = 3t and different m a . In this case 
£ = 1. DOS was calculated using Eq. (Iff [i where the integral 
over p £ [—1, 5] was discretized using 1000 points. Notice the 
different scale (right) for the inverse localization lengths. 



transmission band in the x direction for E G [— 2t, 2t] as 
long as g(p — 0)^0. 

The matrix P has the structure of a ID Anderson 
model with disordered onsite potentials, j3 y , constant 
hopping strength, £ y = £, and hardwall boundaries. By 
Gershgorin's circle theorem. ) 61 ' 62 the spectrum of eigen- 
values p as L y — > oo fills the interval p e [— 2|£| + mp — 
Wp/2,mp + Wp/2 + 2\^\\. The condition for the existence 
of a perfect conducting channel is then 



iei> 



4 ' 



(22) 



The above inequality can be satisfied either by changing 
£ or the /3 y distribution. As an example, we see that when 
mp = 0, the resonant channel will emerge (for large L y ) 
for any value of £. Numerical results for DOS and local- 
ization length for xy-systems in the thermodynamic limit 
(L x = L y = oo) when the above inequality is fulfilled are 
shown in Fig. [B^b). We emphasize that for large L y , 
due to the nature of the DOS in the ID Anderson model 
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FIG. 7: (Color online) Evolution of det P (right scale) and lo- 
calization length X X (E — 0) (left scale, black solid and dashed 
lines) as a function of £ for the a;i/-system used in Fig. [S] Four 
resonant values appear in the interval shown £ = 0.26042, 
0.41092, 0.82164, and 1.2992. The shaded regions highlight 
the intervals of £ corresponding to ±5% and ±10% deviation 
around the resonances. The dashed line corresponds to cal- 
culations of X X (E — 0) where each y XyV value was allowed 
to deviate randomly up to ±5% from the target 7^ = 
value. 



(smooth and differentiable), the condition g(p = 0) 7^ 
implies that the fraction of values of p arbitrarily close 
to is finite, and therefore the number of channels with 
very large localization lengths grows with L y . These will 
then lead to a dense continuum of effectively extended 
states for E G [— 2t, 2t] in systems with large but finite 
L x and L y whenever the above condition is satisfied. 

We can estimate the relation between X x and L y for 
nearly resonant channels, i.e. when (|22[) holds. If we as- 
sume that g(p) = g is roughly constant around p = 0, 
then for finite L y we will find a value of p as small as 
\p s \ ^ {2L y g)~ l . The localization in the corresponding 
decoupled channel — which provides the maximum local- 
ization length — , characterized by a disorder distribution 
of width |p s |W a and mean p s m a , at E = is given by 
(cp. Fig. E| 



A*(.E = 0) 



24 



[At 2 



since (pslVFa is smalL 57 ' 63 This leads to 



2 2 1 



384L 2 u(? 2 

K{E = 0) > 



24m 2 



(w a /ty wi 



(23) 



(24) 



Therefore the minimum localization length at the band 
centre of the nearly resonant channels, which will emerge 
when Eq. (|22j) is fulfilled, scales as L 2 . We can roughly 
estimate the order of magnitude of g by assuming a con- 
stant DOS which leads to g = (Wp + 4|^|)' 1 . For the 
typical parameters used in our simulations, we then ob- 
tain X X (E = 0) > 3L 2 as approximate lower bound. We 
note that this estimate for large L y already works quite 
well for the L y = 10 case of Fig. [7] 
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D. Effective delocalization in y 

In general, for large enough L yi localization of the 
wavefunctions given in (|13[) is to be expected in the y- 
direction, with localization lengths determined by the 
eigenvectors of P (cp. Fig. [TJ). The properties of the 
ID disordered system described by the latter matrix de- 
pend on £ and /3 y , corresponding, respectively, to hopping 
and on-site energies. Hence we can write the localization 
length in the y-direction as 



24 



4£ 2 - (p - mp\ 



(25) 



in the weak disorder approximation, |£| ^> Wp/^/l2 - 57]63 
As explained in Section Hi El the transverse localization is 
determined by the value of p, and thus all states of a cer- 
tain decoupled channel characterized by p c [see Eqs. ©] 
will exhibit the same y-spreading. Localization in y can 
then also be tailored by changing £ and/or the disorder 
distribution /3 y . For example, for fixed disorder, X y can 
be increased by choosing larger £ values. 

Expression (|2"5j) can be used to estimate the value of 
£ that we need if we want the extended states that we 
have generated in the rr-direction to be also effectively 
delocalized (\ y > L y ) in the vertical direction. By using 
Eq. (f2~5j) for p = we estimate that these states will be 
effectively delocalized in the y direction when 



1/2 



(26) 



For the case rap = we see that for systems a few hun- 
dred sites wide a value of £ ~ Wp will roughly ensure 
effective delocalization over the whole system of the reso- 
nant channel states, as shown in Fig.[5Ja). Alternatively, 
we can also find effective delocalization in y, while keep- 
ing the localized character in the x direction, as displayed 
in Fig.HJb). 

Therefore, the localization properties of xy-disordcr 
can be engineered almost at will, and independently in 
the x and y directions, giving rise to the 4 basic configu- 
rations shown in Figs. Q] and [5] 



IV. ENGINEERING OF TRANSPORT AND 
SPECTRAL PROPERTIES 

A. Transport regimes in xy-disorder 



from transmission probabilities in x and y with open 
boundary conditions — averaged over 100 realizations of 
the system, displayed by the color density plots in Fig. |H1 

Transmission in the x direction emerges from the exis- 
tence of resonant or quasi-resonant channels, whereas in 
y it is due to effectively delocalized states only. Thus it 
must be clear that in the thermodynamic limit, L x , L y — > 
oo, only the efficient transport along x remains. The dia- 
gram is system-size dependent with respect to the trans- 
port properties in the y-direction. Nevertheless, for finite 
systems, there exists the interesting possibility of realiz- 
ing switching devices whose transmitting behaviour in x 
and y is highly controllable. 

Different configurations of the diagram can be obtained 
by changing the disorder parameters W a , m a , Wp and 
rap. For example, rap controls the appearance of the res- 
onant channels in x, as well as the width of the plateau 
region in the border of effective delocalization in y. On 
the other hand, ra a ^ breaks the symmetry of the dia- 
gram, the region of efficient transport in y shifts laterally 
preserving the plateau region and decaying asymmetri- 
cally on the sides. In Fig. [TOl we show an alternative 
diagram of mp vs E, where different transport regimes 
can be triggered simply by changing the mean value of 
the j3 y distribution, even for a fixed disorder realization. 



B. Coexistence of extended and localized states 

In a finite a;y-system, the fulfilment of the resonance 
condition (as discussed in Scction lHI Ajl induces L x Bloch 
states in the x direction for E 6 [— 2t, 2t\. Nevertheless, 
the other (L y — 1)L X states coming from the remaining 
decoupled channels will be localized in x, and can in prin- 
ciple have energies inside the range [— 2t, 2i\. Therefore, 
in general, we should expect coexistence of extended and 
localized states in [—2t, 2t\. The overlapping of the spec- 
trum of resonant and localized channels can be seen in 
Fig. [H^a). We have checked that the level spacing dis- 
tribution in this region includes that for Bloch states on 
a background of Poissonians corresponding to localized 
states, as also found in Ref. for a ladder model sup- 
porting coexistence. 

As L x and L y grow, the range [—2t, 2t] will poten- 
tially become densely populated by extended and local- 
ized states. The filling of that energy interval by Bloch 
states in a continuous manner will render the presence 
of localized states irrelevant, as far as transport in the x 
direction is concerned. 



In Fig. [9] we show the diagram of transport regimes 
|£| versus E for a 2D system of size L x = L y = 1000. 
By comparing the magnitude of the localization length 
against the system size, we can distinguish regions of ef- 
ficient transport in direction x, in y or in both. The 
lines separating the different phases are calculated in Ap- 
pendix [Cj We confirm the validity of the diagram by nu- 
merical calculations of localization lengths — estimated 



C. Avoiding the coexistence of extended and 
localized states 

The c-th decoupled channel is characterized by an 
on-site energy distribution with mean p c m a and width 
|p c |Wa, and whose spectral boundaries are given by 
Eq. (fTo| . Thus m a controls the position of the spec- 
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FIG. 8: (Color online) Probability density l"!/^,^ 2 of eigenstates of an ^-disordered system effectively delocalized in y with (a) 
Aa; > L x and (b) \ x < L x . The parameters characterizing the potential are m a = 0, W a — 3t, mp = 1, Wg = 2, £ y = £ = 2 and 
eigenenergies (a) £?/t = 1-7535, (b) -E/i = —0.008711. The states belong respectively to decoupled channels with (a) p — 0.044 
and (b) p = —1.2029. Notice that £ is chosen to satisfy Eqs. (|22|) and (|26|) only. In both cases the system length and width 
are L x = 100 and L y — 50, respectively. The left and bottom panels show, respectively, the distributions |xy| 2 an d \4>x^ that 
construct the eigenstate according to Q13p . The color is determined by the value of — l°g(L) IV , ».v| 2 where (L) = L x L y (in 
this scale, 2 is the average value of an extended state). A second order interpolation of the distributions is used to smooth the 
visualization. 



trum of the different channels. The avoidance of coexis- 
tence is only possible if the a x 's have all the same sign, 
i.e. \m a \ > W a /2. Otherwise the spectrum of all de- 
coupled channels always includes E = 0, as the lower 
and upper bounds have opposite signs, thus giving rise 
to overlapping with the range [— 2t, 2t]. 

For simplicity wc consider m a > 0. Then we can work 
with | p c | and consider positive energies only without loss 
of generality. In order to avoid coexistence, we have to 
ensure that the lowest energy of all regions of localized 
states is greater than 2t. In view of Eq. (|16|) this condi- 
tion is written as 

(m a -^jmin\p c \ > 4t, (27) 

where min|p c | stands for the minimum non- vanishing 
\p c \. Let us recall that the resonance condition is sat- 
isfied and thus zero is an eigenvalue of P [Eq. (JBJ]. 

Once the disorder realization of the system is gener- 
ated, i.e. fixed sequences {cv^}, {(3 y } and £, we can ob- 
tain the value of min \p c \ and then change m a to satisfy 
the inequality. We note that ([27]) is formally independent 
of the system size, but only reasonable for moderate L y . 
Since we are in the region of resonance, for large L y the 
min \p c \ will get arbitrarily small (Section IIII C|) and the 
required m a to avoid overlapping will be very large. 

It is also possible to obtain some general bounds for 
m a without having to apply the above inequality to each 
particular realisation of the disorder. For a fixed value 
of the width of the system L y , we can estimate min \p c \ 
from the distribution g(p) around p ~ 0. A rough es- 
timation is obtained by assuming a constant distribu- 
tion, g = (Wp + 4|£|) -1 , from which it follows that 



min|p c | ^ (W/3 + 4|£|)/L r Condition $27} translates 
into 

which gives an approximation, usually overestimated, of 
the minimum m a value to avoid coexistence. In Fig. [TT] 
we show examples of spectrum engineering of the over- 
lapping between extended and localized states. 

V. CONCLUSIONS 

Our results show how a controlled choice of quite a 
simple, structured disorder as given by Eqs. (j4]) and ([5]) 
can lead to extended states. By varying the parameters 
of the individual disorder distributions, we can engineer 
the localization and transport properties of the system 
independently in each direction. We can furthermore 
control the energy coexistence of extended and localized 
states in the spectrum. Although our approach is based 
on the Anderson model, we expect our results to hold 
in the many realms of Anderson localization physics in 
g eneral»M-70 

For an experimental realisation of our model, the ro- 
bustness of our results to small deviations away from 
the exact resonance conditions for quasi-lD systems, as 
shown in Section MI B[ and the broadening of the res- 
onance condition for 2D systems (Section IIII CP is of 
course reassuring. Nevertheless, the disorder does re- 
quire single-site control in order to adjust the e XtV and 
values as needed. The key feature expected of any 
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FIG. 9: (Color online) Diagram of transport regimes for |£| 
versus £ in a 2D system with xy-disorder of size L x — L y — 
1000 and parameters W a = t, m a — 0, Wp = 1, mp = 2.5. 
The dashed lines mark the limits of the efficient transport 
regimes obtained analytically and given by Eq. (|22[) (straight 
horizontal dashed line), Eq. (|C2[I (vertical dashed lines) and 
Eqs. (|C6[) and (|C7b[) (bent dashed line). The density plots in 
the background (color scales are indicated above the main 
panel) show numerically calculated localization lengths \ x 
(white-to-blue semitransparent) and X y (white-to-red). The 
(E, |£|) space was discretized in 675 points (27 x 25) for the 
calculation, and the results averaged over 100 disorder real- 
izations at every point. The density plot uses a linear interpo- 
lation to smooth the visualization. Labels x and y highlight 
regions where > L x and \ y > L v respectively. The gray 
wedges at the lower corners correspond to regions outside the 
spectrum of the system. 
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FIG. 10: (Color online) Diagram of transport regimes for mp 
versus E in a 2D system with xy-disorder of size L x = L y — 
1000 and parameters W a = t, m a = 0, Wp = 1, £ = 3.4. The 
dashed lines mark the limits of the efficient transport regimes 
and can be obtained from the equations in Appendix [C] The 
density plots in the background show numerical calculations 
of the localization lengths as in Fig. [9] 
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FIG. 11: (Color online) Spectrum of zy-system with resonant 
channel and parameters L y — 10, W a = 2t, £ = 4.23152 
for different m a . Black lines indicate DOS for L x = oo as 
obtained from FEF, while blue (red) points on the x axis 
mark eigenenergies of extended (localized) states for a finite 
realization with L x = 100. The shaded gray region highlights 
the contribution of the resonant channel to the total DOS. 
The inset shows the /3 y sequence considered characterized by 
Wp ~ 3. According to Eq. ((281) blue-red coexistence should 
be avoided for m a > 3. 



experimental endeavour to measure our model is hence- 
forth the attempt to achieve full control by implementing 
single site resolution for quantum state manipulation as 
well as quantum state analysis. Fortunately, important 
progress in this direction has already been achieved, e.g. 
in optical lattices42r— A scheme combining a disordered 
optical potential and an additional harmonic trap has 
been recently proposed^ to observe experimentally the 
coexistence of extended and localized wavefunctions. We 
emphasize that our xy-disordcr can provide a genuine and 
controllable coexistence of extended and localized wave- 
functions that does not occur in usual disordered systems 
exhibiting dclocalization transitions. 

While we have analysed 2D models for simplicity, the 
product structure for the disorder presented here can be 
straightforwardly extended to 3D systems, where similar 
features are to be expected, in particular, the possibility 
to engineer the transport properties in the three spatial 
directions independently 
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Equation (|B1[) can be expressed in standard transfer- 
matrix form as 



t-\El-e x ) -1 
1 



(B2) 



Appendix A: Functional equation formalism 

The DOS and the localization length of an infinite ID 
disordered chain, can be calculated numerically using the 
FEFi££ This technique has been applied to obtain the 
spectral properties in the thermodynamic limit of differ- 
ent quantum-wire models with correlated and uncorre- 
cted disorder.— ~— Here we give a brief overview of the 
formalism applied to ID tight-binding models with diag- 
onal disorder, described by the equation (E — e x )<f) x = 
4> x +i — <f> x -i) where the on-site energies are obtained ran- 
domly from a continuous distribution V(e). For such a 
system, the functional equation reads, 

W{9) = J deV(e)W(T(9; E,e)) - VV(tt/2) + 1, (Al) 

where T(6;E,e) — arctan (E — e — 1/ tan#), with the 
conditions 9 € [0, n) and W{9 + nir) = W(9) + n for inte- 
ger n. The angular variable, 9, follows from the represen- 
tation of the wavefunction amplitude <j> in polar coordi- 
nates. The function W(9) corresponds to the cumulative 
distribution function of 9 in the disordered chain in the 
thermodynamic limit. The DOS g(E) and the Lyapunov 
exponent T] X (E) as functions of the energy are written in 
terms of W(9) as, 



9(E) 



dW(n/2) 



dE 



(A2) 



Vx(E) = l - J deV(e)\nf(ir;E,e) 



d£r{£) f mE, £ e) > (A3) 



where f(9,e) = 1 - (E - e) sin26» + (E — e) 2 cos 2 9, and 
the prime indicates differentiation with respect to 0. By 
solving numerically Eq. (|A1[) for different values of the 
energy, the spectral and localization properties of the in- 
finite system can be obtained using the latter expressions. 



Appendix B: Transfer-matrix method 

The transfer-matrix method (TMM) allows for a very 
memory efficient way to iteratively calculate the decay 
length of electronic states in a quasi-lD system of width 
L y for lengths L x 3> L y J^ Equation ^ has to be rear- 
ranged into a form where the ^ x +i amplitudes of sites 
in layer x + 1 — when x is chosen as the direction of 
transfer — is calculated solely from parameters of sites 
in previous layers x and x — 1, 



where and 1 are the L y x L y zero and unit matrices, 
respectively. Formally, the transfer matrix T^. is used to 
'transfer' electronic amplitudes from one slice to the 
next and repeated multiplication of this gives the global 



transfer matrix tl x = 



T, 



r = lini£ 

ues e ±v \ 



-t 



The limiting matrix 

t Lx t l J exist a 76 i 77 and has eigenval- 

i = 1, . . . , L y . The inverse of these Lyapunov 
exponents iji are estimates of decay/localization lengths 
and the physically most relevant largest decay length in 
the x direction is X x = l/minrjj. 

Appendix C: Transport regimes 

The transport properties of the system along x and y 
can be engineered by tuning the parameters of the sys- 
tem: m a , m/3, Wq,, Wyg, £ and E. A qualitative diagram 
of the transport regimes, for sufficiently large but finite 
L x and L y , can be obtained by using the weak disorder 
expressions for the localization lengths and X y . For 
simplicity we assume here that m ai mp ^ and we only 
focus on relations |£| versus E. 

a. Transport along the x-direction: The basic con- 
dition to have resonant channels for large L y is given by 
Eq. (|22p . Let us approximate that g(p) = g is constant 
for all p and thus g = (Wp + 4|£|) -1 . Then for finite 
L y , when the above condition is satisfied, we will find 
a value of p as small as \p s \ < (Wp + 4\£\)/(2L y ). The 
maximum localization length in x (provided by a chan- 
nel with a disorder distribution of width |p s |W Q -C 1 and 
mean p s m a ~ 0), then satisfies 



96Ll(At 2 -E 2 ) 



E G [-2t,2t] 



(CI) 




x+l 



t-\Et-e x )* x -* a 



(Bl) 



w^wp+mr' 

We can ensure efficient transport in x by enforcing that 
the minimum bound for X x is larger than L x . This im- 
plies 



(C2) 



This condition gives a continuum region of £ values pro- 
viding transmitting behaviour in x. Note, however, that 
X x > L x will also be achieved outside this region, around 
resonant values of £ (p s =0). 

b. Transport along the y-direction: The localization 
length in y is determined by the localization properties of 
the eigenstates of the matrix P. The localization length 
in p is given by 

24 r i 
X v(P) = TI72 ^ 2 -(p-mpf , (C3) 
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for p e [-2|£| + m p ,mp + 2|£j] and |f| > Wp/y/tt. In 
this weak disorder approximation the localization length 
has a parabolic dependence in p with its maximum value 
at p = mp. Since the localization length in y is deter- 
mined by the value of p, all states of a certain decoupled 
channel will exhibit the same y-spreading [see Eq. fl5J]. 
The spectrum of a channel characterized by p c is 



E € 



-2t 



crW a 



,2t+p c 



(C4) 

where cr = sign(p c ). For a given energy £7, the maxi- 
mum localization length in y, and thus the relevant for 
transport processes, is given by X y (p c ) such that p c is the 
closest value to mp whose spectrum includes E. 

For p c = mp we get the maximum attainable value of 
\ y in the system, and thus X y (E) = 96£ 2 /Wj for 



(E + 2t) I (m a — aW a /2). The expression and its energy- 
range of validity depends on whether the contributing p c 
is positive (cr = 1) or negative (cr = —1), as well on the 
sign of (m a — W a /2). 

Therefore, from the plateau structure in the interval 
(|C5j) . the localization length in y decays with E. For 
m a = the dependence of X V {E) is symmetric around 
E = 0. A non- vanishing m a induces a shift in the posi- 
tion of the plateau and can lead to an asymmetric decay 
of \ y {E) from its maximum value, left and right from the 
plateau. These features can be observed in the numerical 
results shown in Fig- |3](b) . 

In order to have efficient transport along y we require 
X y Ly. Following the reasoning given above it is pos- 
sible to write general relations for |£| to ensure the latter 
condition. In order to avoid cumbersome expressions, we 
show them here only for the case m a = 0: 



E 



-2t 



mp 



2 



, 2t + mp 



Wo 

2 



(C5) 



in Eq. (|C3[) to obtain X y as function of the energy 



For energies smaller than the interval (|C5[) . it can be 
seen that the maximum X y is attained for p c such that E 
coincides with the lower edge of the spectrum of the p c - 
channel. Thus X y (E) is obtained from Eq. (|C3[) for p c = 



|f|> 



4V6 ' 



E e 



-2t - mp — ,2t + mp — 



since any other p c whose spectrum overlaps with this an( j 
range will provide smaller values of X y . The localization 
length in y thus displays a characteristic plateau struc- 
ture [see Fig. [3(b)] whose width is proportional to the 
parameter mp. 

For energies larger than the interval above, there is 
overlapping of the spectra of channels with p > mp (recall for 
that we assume m a , mp 0). For a given energy E, the 
maximum X y is attained for p c such that E coincides with 
the upper edge of the spectrum of the p c -channel, since 
this will be the value of p closest to mp. Therefore we 
can identify p c = (E — 2t)/(m a + W a /2) and substitute 



' LyW$ 

24 



\E\ - 2t 
W a /2 



mp 



(C6) 



(C7a) 



\E\e 



2t + mp^,2t+(mp + m)^- 



(C7b) 



The latter relations arc symmetric around E — 0. As dis- 
cussed above, a non-vanishing m a will shift the plateau 
structure (|C6[) . and it can break the symmetric behaviour 
around it. The formulas above provide a qualitative un- 
derstanding of the influence of the different parameters 
on the transport regimes in the y direction. 
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